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Abstract 

Exact analytic solution for the probability distribution function of the non-inertial 
rotational diffusion equation, i.e., of the Smoluchowski one, in a symmetric Maier- 
Saupe uniaxial potential of mean torque is obtained via the confluent Heun’s func¬ 
tion. Both the ordinary Maier-Saupe potential and the double-well one with vari¬ 
able barrier width are considered. Thus, the present article substantially extends 
the scope of the potentials amenable to the treatment by reducing Smoluchowski 
equation to the confluent Heun’s one. The solution is uniformly valid for any barrier 
height. We use it for the calculation of the mean first passage time. Also the higher 
eigenvalues for the relaxation decay modes in the case of ordinary Maier-Saupe 
potential are calculated. The results obtained are in full agreement with those of 
the approach developed by Coffey, Kalmykov, Dejardin and their coauthors in the 
whole range of barrier heights. 
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1 Introduction 


Rotational reorientations take place in liquid crystals, polymers, proteins, 
lipids and many other systems. They are investigated by means of variety of 
experimental techniques such as NMR, dielectric relaxation spectroscopy, flu¬ 
orescence depolarization, etc (see, e.g., [1] and refs, therein). The basic tool for 
their theoretical description is the Fokker-Planck equation [2]-[5]. For the case 
when the inertial effects are negligible (overdamped limit) the latter is reduced 
to the Smoluchowski equation (SE). There is vast literature on their study and 
applications [2]-[46]. The applications of them to the dielectric spectroscopy 
[29], [30], [26], fluorescence depolarization [16] and NMR relaxation of liquid 
crystals (see [31], [32], [33] and refs, therein) have been thoroughly explored. 
Also, the extension of the theory from ordinary diffusion to the fractional one 
is intensively studied (see [25], and refs, therein). Besides SE arises in the 
theory of ferromagnetism (where it is called Brown’s equation) [34], [35], [44]. 
Both in the theory of molecular rotational motion in a uniaxial potential [45], 
[27], [25] (occurring for dynamics of liquid crystals, dielectric relaxation, etc.) 
and in the theory of ferromagnetism [44] the so called Maier-Saupe potential 
is widely used and by right considered to be a basic model. 


The most powerful approach to the problem was developed by Coffey, Kalmykov, 
Dejardin and their coauthors in the above mentioned papers [7], [17], [25], [26], 
[27], [28]. For the sake of brevity we further call it as CKD. The approach is 
based on the expansion of the probability distribution function as a series of 
spherical harmonics. This method is well suited for the potentials of mean 
torque that can be expanded in terms of spherical harmonics. It results in 
an inhnite hierarchy of differential-recurrence relations for the moments (the 
expectation values of the spherical harmonics). On the one hand CKD is an 
exact approach that uses no approximations and imposes no physical limi¬ 
tations. The undoubted merit of CKD is that it is very general and can be 
applied to enormous variety of different physical situations described by cor¬ 
responding types of potentials. On the other hand this approach makes use 
of foreign eigenfunctions (Legendre polynomials) for SE under consideration. 
For this reason, a complementary approach was suggested by the author [43] 
that yields an exact solution in purely mathematical sense. The latter uses 
the expansion of the solution over own eigenfunctions of the Smoluchowski 
operator rather than that over Legendre polynomials. The own eigenfunctions 
are expressed via the confluent Heun’s function (CHF). The advantage of 
our approach is that it enables one to avoid dealing with inhnite hierarchies 
of differential-recurrence relations that is a complex and laborious problem 
in itself. Besides it yields directly the probability distribution function itself 
rather than statistical moments (i.e., values averaged over it). Such option is 
convenient for, e.g., applications of the theory to relaxation processes in NMR. 
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However the approach of [43] was developed for a rather specihc potential 
(double-well Maier-Saupe one with fixed barrier width parameter). The main 
aim of the present paper is to show that our approach has wider applicability. 
We substantially extend the range of potentials amenable to treatment by this 
method and apply it to both the ordinary and the double-well Maier-Saupe 
potentials while the latter may have variable barrier width parameter. The 
consideration is carried out within the framework of a uniform approach based 
on the general Maier-Saupe potential of which both ordinary and double-well 
options are merely particular cases. Nevertheless, despite the above mentioned 
extension of the options our heuristic potential is still a symmetric uniaxial 
Maier-Saupe mean field one 

= a sin^ — 2c In sin 

ksT 


where a is the barrier height parameter and c is the barrier width parame¬ 
ter that takes the values from the range 0 < c < cxd. The variable barrier 
width parameter c enables one to eliminate the logarithmic term at wish. The 
standard uniaxial potential 


ksT 


a sin^ ip 


is merely a particular case c = 0 of the previous one so that both cases are 
amenable to exact analytic treatment via the CHF in a uniform manner. The 
particular case c = 0 was explored in details by CKD in [45], [27], [25]. The 
results obtained there are the reliable reference data for the present approach. 
The SE for the standard form of the potential c = 0 can be as well solved 
exactly in terms of known functions in the same way as that for the potential 
with the logarithmic term c > 0. From the purely mathematical point of 
view it seems natural to treat both cases simultaneously. The particular case 
with fixed barrier width parameter c = 1 explored in [43] is distinguished by 
extreme simplicity of the boundary conditions for the characteristic equation 
for the spectrum of eigenvalues. 

From the physical point of view the potential with the logarithmic term is con¬ 
sidered for the following reason. For rigid rotators without internal motions (a 
nematic liquid crystal molecule or a hne single-domain ferromagnetic particle) 
the minima of a uniaxial potential over the polar angle ip in the laboratory 
frame are at '^ = 0 and = tt. In these cases the standard Maier-Saupe poten¬ 
tial is relevant and there are no reasons to expect that the mean held theory 
can provide the appearance of the logarithmic term in such systems. In con¬ 
trast the latter hnds visible physical motivation within the context of internal 
(segmental) motions in polymers and proteins (e.g., those of hi-loops). The 
minima of a uniaxial potential over the polar angle tp in the molecular frame 
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for a segment of structure are not necessarily at '^ = 0 and if) = n and we need 
a genuinely double-well potential to describe such motions. The logarithmic 
term in the potential is a simplest way to shift the equilibrium positions from 
= 0 and ifj = n. However an important difference should be stressed. The 
segmental motions can not proceed with the isotropy of the potential over the 
azimuthal angle. In most cases they proceed with a hxed azimuthal angle that 
requires a different normalization of the probability distribution function than 
for the case of isotropy (see [43] for more details). Otherwise the considera¬ 
tion follows a similar line as for the case of isotropy. In the present article we 
choose the normalization with the isotropy of the potential over the azimuthal 
angle for the sole reason: to provide the possibility of direct comparison of our 
results with those of the approach developed by CKD. We show that our ap¬ 
proach is valid for this case and hence can be sure that it will remain valid for 
a realistic case of the potential with a hxed azimuthal angle. 

As mentioned above our approach yields directly the probability distribution 
function for SE. The latter enables one to obtain any required average in a 
uniform manner. We calculate with its help the mean hrst passage time to 
demonstrate that our approach leads to identical results with those of CKD. 
We pay special attention to the case of ordinary Maier-Saupe potential be¬ 
cause for it there is a full set of data from CKD. In our opinion it is of interest 
to compare the results of the present approach based on reducing SE to con- 
huent Heun’s equation with those of CKD for this case. On the one hand 
as mentioned above the ordinary Maier-Saupe potential is a basic theoretical 
model for materials based on both nematic liquid crystals and ferromagnetics. 
On the other hand the case of ordinary uniaxial Maier-Saupe potential de¬ 
spite of its simplicity is the cornerstone basic example for the theory. Hence, 
any discrepancies in this key issue are impermissible. We show that there is 
no divergence between the results of two approaches in this case. There is an 
important issue where our approach can be especially helpful. In [41], [42] the 
problem of higher eigenvalues for the relaxation decay modes was raised. It is 
related to the question whether the relaxation proceeds via a single mechanism 
or not. Within the framework of CKD the calculation of higher eigenvalues 
requires laborious numerical work by the method developed in [42]. In our 
approach the eigenvalues are the roots of the a characteristic equation for the 
corresponding boundary problem for the Smoluchowski operator. This equa¬ 
tion includes CHE and its derivative. Maple makes Ending its roots a routine 
procedure in all practically important cases. However, some technical problems 
arise at calculating higher eigenvalues (e.g., A 5 ) by Maple for small barriers. 
They presumably result from known deficiencies of realization of CHE deriva¬ 
tive in this symbolic computational software package. We discuss the issue 
below. 

The paper is organized as follows. In Sec. 2 the problem under study is for¬ 
mulated. In Sec. 3 the solution of SE is presented. In Sec. 4 the spectrum of 
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eigenvalues is derived. In Sec. 5 the probability distribution function is ob¬ 
tained. In Sec .6 the MFPT is calculated with its help. In Sec. 7 the results are 
discussed and the conclusions are summarized. 


2 Smoluchowski equation 


For a uniaxial potential of mean torque SE under consideration is [45], 

|27], [25] 




_ ^ 

sin Ip dip 


sinipfiip^t) 


dVjip) 

dip 


+ 


1 <9 [ . ,df{ip,t) 

- 7 — sm Ip —-- 

sin Ip dip [ dip 


( 1 ) 


where r is the characteristic relaxation time for isotropic non-inertial rota¬ 
tional diffusion (e.g., td the Debye one for the theory of dielectric relaxation 
or tm the Neel one for the theory of ferromagnetism) and /9 = 1/ {ksT). Its so¬ 
lution is the probability distribution function f{ip,t) that must be normalized 
so that its integral over the whole space gives 


1 

Att 


2tt 


j dip J dip sin ^7 f{ip,t) = - J dip sinip f{ip,t) = 1 


( 2 ) 


The normalized initial condition takes the form 

/(^, 0 ) = —^ S{ip-ipo) 
sm ipo 

We introduce a new variable 
X = cosip 


( 3 ) 


( 4 ) 


Then the equation takes the form 

df{x,t) 1 d 

dt 2 t dx 

where the dash means the derivative over variable x. 


(1 


X 


df{x,t) 

dx 


+ l3f{x,t)V\x) 


( 5 ) 
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Further we consider a general Maier-Saupe uniaxial potential of mean torque 
potential 


r ( i ) = ^ ^ '‘*(1 “ 


( 6 ) 


The ordinary Maier-Saupe uniaxial potential of mean torque potential consid¬ 
ered in [45], [27], [25] corresponds to the particular case c = 0. The double-well 
Maier-Saupe potential with fixed barrier width parameter suggested in [43] 
corresponds to the particular case c = 1. Both barrier height E and width w 
depend on a and c 



Inversely one obtains 


a = I3E 




-1 


w 


c = ^E\l-- 



( 8 ) 


For the ordinary Maier-Saupe potential (c = 0) we have w = 2. The de¬ 
pendence of the potential (6) on the parameters is depicted in Fig.l. The 
requirement for the potential to be a double-well one is c < a. 


3 Solution of Smoluchowski equation 


Equation (5) can be solved by separation of variables 
= I 7 (x)exp (-pt) 


(9) 


The equation for the function r]{x) is 

(l - 7]’’^{x) - 2x [a (l - x^) - (c - 1)] v'xix) + 


pr + {c — a) + 3ax‘^ rj{x) = 0 


( 10 ) 
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y{mw) 




tp (radians) 


Fig. 1. The symmetric Maier-Saupe uniaxial potential of mean torque 
V{ip)/ {ksT) = a sin^ V’ “ c In {sit? if). For c > 0 it is called the double-well 
Maier-Saupe potential to distinguish it to the ordinary Maier-Saupe one c = 0. The 
value of the parameter u is ct = 14. 

We introduce a new variable 

y = x‘^ ( 11 ) 


The equation takes the form 


y{y - l)^ra(2/) + 


-^y"^ + o (1 + 2cr - 2(c - 1)) 1 / - - 


iviy)- 


- [pr + (c - a) + ?>ay] r]{y) = 0 


( 12 ) 


It belongs to a class of the so-called confluent Heun’s equation [47]. Equa¬ 
tion (12) has fundamental solutions that can be expressed via the confluent 
Heun’s function (CHF). The latter is a known special function [47], [48], [49]. 
At present it is realized explicitly in the only symbolic computational soft¬ 
ware package Maple as HeunC and its derivative HeunCPrime (see [48] for 
detailed discussion of the merits and drawbacks of this computational tool in 
Maple). Returning to the variable x we can write the general solution of (10) 
as follows 


ri{x) = Crj^^fx) + Dp^'^fx) 


( 13 ) 
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where 


ri^^\x) = HeunC^^'> -c, (c + ^ , 

7]^^\x) = xHeunC^'^^ |^-a, -c, (c + 


(14) 

^;^^)(i5) 


4 Spectrum of eigenvalues 


The boundary conditions are obtained from (10) by setting x = ±1 


2x{c — l)ri'^{x) + 2 [pr + c + 2a] r]{x) 



±1 


(16) 


The characteristic equation for our boundary problem is split apart to a pair 
of equations (omitting for the sake of brevity the parameters of CHFs and 
their derivatives) 

2(c — l)HeunCPrime^^^ (; 1) + (pr + c + 2a)HeunC^^'^ (; 1) = 0 (17) 

2 (c — VjHeunCPrime^"^^ (; 1) + (pr + 2c — 1 + 2a)HeunC^‘^^ (; 1) = 0(18) 

In case c = 1 they are reduced to simpler ones [43] 


HeunC^^'^ (; 1) = 0 (19) 

HeunC^'^'^ (; 1 ) = 0 ( 20 ) 

We denote 

= P^n^r ( 21 ) 

llin = (22) 


The equations (17) and (18) can be solved only numerically but Maple copes 
with this problem. One always obtains Ai = 0 while the lowest Oi is always 
nonzero and the corresponding term mainly determines the behavior of the 
distribution function f{'ip, t) and MFPT calculated with its help. The behavior 
of Oi is depicted in Fig.2. Higher eigenvalues are plotted in Fig.4 and Fig.5. 



In(ni) 



|-#-c=5 present paper approach 


c=3 present paper approach 


|^c=l present paper approach 


IB c=0 present paper approach | 
|— Coffey et al. approach for ln(AiT)| 


Fig. 2. The dependence of the first non-zero eigenvalue for the boundary problem of 
rotational Smoluchowski equation on the parameter a. The latter dehnes to a greater 
extent the barrier height of the double-well Maier-Saupe potential. The continuous 
line is the result of the approach developed by Coffey, Kalmykov, Dejardin and their 
coauthors (32) for the particular case of ordinary Maier-Saupe potential. 

5 Probability distribution function 


We return to the variable '0 and write the solution of SE (1) as 
fi'ip, t) = Cir][^^ (V') + Di exp (-fiit/r) ^ (V')+ 

C2exp (-A2t/r) + D2exp (-122^/^) + ... (23) 


The coefficients Cn and Dm are calculated with taking into account the fact 
that being the solution of the boundary problems and are 

orthogonal of each other 

TT 

J dxl) sinV' Vn\'4’)Vm{'4’) = ^mnSkl (24) 

0 


where k,l = 1,2 and n, m = 0,1, 2,... . At time f = 0 (23) yields (with taking 
into account (3)) 

= Cir][^\'ijj) + Dir]^^\'ijj) -k C2V2 \'ip) + D2r]i^\ilj) + ... (25) 
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Then the coefficients Cn and Dm are obtained by multiplying (25) by 
or 77 )^) ('?/’) and integration over tjj with taking into account (24) 




2hnH^0) 


TT 

J dtlj simjj 


LO 



and 


D 


ra 


2hmH^o) 


TT 

J d'lp simjj 


LO 



(26) 


(27) 


Thus, we have the explicit algorithm for obtaining the coefficients Cn and Dm 
along with the corresponding eigenvalues and Dm that makes the solution 
of our problem to be completed. Substitution of all these values into (23) yields 
the required probability distribution function. 


6 Mean first passage time 


To exhibit how the result obtained may be useful we calculate with its help 
the dependence of MFPT on the parameters of the potential a and c. Our 
effective Maier-Saupe uniaxial potential of mean torque is 


1 


V{ijj) = a sin^ ijj — c In (^sin^ il?j 


(28) 


It has two local minima the left of which we denote as i/jl and the maximum 
at ipM = 7 J'/ 2 . Following Risken [ 6 ] one can obtain for the non-averaged MFPT 
Ti that characterizes the ability for the system to reach the barrier top ipM 
starting from the left minimum jpQ = the expression (see, e.g., [43] for 
details) 


Ti(^o) 


dM 


dtjj simjj dT T 


dT 


(29) 


After straightforward calculations we obtain 


bM „ 

TijpJ^) = J dpj sin 77 |^h®(^) + + •••j (30) 
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ln«ri>/T) 



0 c=5 present paper approach 
^ c=3 present paper approach 
^ c=l present paper approach 
B c=:0 present paper approach 

— Coffey et al. approach I 
Coffey et al. approach II 

— Coffey et al. approach II 
Coffey et al. approach II 


Fig. 3. The dependence of the mean first passage time for the transition from the 
well of the symmetric double-well Maier-Saupe uniaxial potential of mean torque 
V{ip)/ [ksT) = a sin^'0 — c In (sin^ on the parameter a. The latter dehnes to 
a greater extent the barrier height of the double-well Maier-Saupe potential. The 
continuous line for the particular case of ordinary Maier-Saupe potential (c = 0) 
is the result of the full-fledged (time dependent) approach of Coffey, Kalmykov, 
Dejardin and their coauthors (33). The continuous lines for cases with c > 0 are the 
result of the restricted approach of Coffey, Kalmykov, Dejardin and their coauthors 
(34) within the framework of the stationary approximation. 


Finally, we calculate the averaged MFPT <Ti > for the case when the initial 
condition is thermodynamically averaged over the left well 


< Ti >= 




sin Tpo exp 



(31) 


It is worthy to recall here that the coefficients Cn and Dm are the function of 
the initial condition 'i/jq. The results of numerical calculation for the averaged 
MFPT < Ti > are presented in Fig.3. 
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7 Results and discussion 


In Fig.2 the first non-zero eigenvalue for the ordinary Maier-Saupe potential 
(c = 0) is compared with the result of CKD for Air given by the empirical 
equation [25], [44] which is valid for all barrier heights 


Air 




2 ^ 3/2 

vT (1 + fr) 


(32) 


Our results are in excellent agreement with those of CKD for the case of 
ordinary Maier-Saupe potential (c = 0). 


In Fig.3 our results for MFPT for the ordinary Maier-Saupe potential (c = 0) 
are compared with those of CKD [25], [44] 


< Ti > 


= 4 


r 


0 z 

p. ^ U 4, /> 


dz 


u 1 — J 

.-1 -1 


dz' e° 


-1 


+ 


[ dz f dz' 


n -1 


.0 


1 - 


-1 


(33) 


We call this formula CKD I. It is obtained within the framework of a full- 
fledged (time dependent) approach taking into account the transient dynamics. 
One can see that our approach is in excellent agreement with CKD I in the 
whole range of barrier heights. The points in Fig. 3 are calculated with the 
only mode Di everywhere but for the values a = 1; 1.5; 2. For the latter we 
have to add the modes D 2 and A 2 to provide agreement with (33). 


There are no results of full-fledged CKD for double-well Maier-Saupe potential 
(c > 0) in the literature. For this reason in Fig. 3 we compare our results for 
this case with those of restricted CKD obtained in the stationary approxima¬ 
tion. The latter in our case of the potential (28) is eq. (1.18.1.7) of [25] 


< Ti > 


r 


= 2 y #' 


a sin^'i/i'—c (n(sin^ i/)') 


siu'^' 


■X 


b' 

J d^J sin^ e-" 

0 

We call it CKD II. In this case the results of two approaches are also in quite 
satisfactory agreement for the whole range where the potential remains to be 
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a double-well one (c < a). Outside this range the very notion of MFPT loses 
its meaning. As earlier all the points in Fig. 3 for the cases c > 0 are calculated 
with the only mode Oi everywhere but for corresponding lowest values of a. 
For the latter we have to add the modes O 2 and A 2 to provide agreement with 
( 34 ). 

Fig. 4 shows that higher eigenvalues are also in excellent agreement with those 
of CKD. The values of A 3 and A 5 for comparison with our results are taken 
from Table. 1 of [42]. However, the implementation of computations within the 
framework of our approach is much more easier than that by the method of 
[42]. It is accomplished with the help of Maple practically automatically. Nev¬ 
ertheless, Maple has problems to hnd the roots for Qm with m > 3 below some 
critical barrier heights, e.g., those for A5 (our Ha) at the values of a = 1 , a = 2 
and a = 3 (see Fig.4). The same phenomenon takes place for with n >3 
as Fig.5 shows. We stress that none of these missed roots is necessary for the 
calculation of MFPT in the agreement with CKD. As was noted above there 
at most H 2 and A 2 are invoked to at some marginally low barrier heights. For 
the particular case c = 1 when the characteristic equation does not include 
HeunCPrime (see (19) and (20)) there was no problem with calculation of 
eigenvalues at any a [43]. Thus, we have drawn the conclusion that the phe¬ 
nomenon may be a result of known dehciencies of realization of HeunCPrime 
in Maple (see [48] for expert opinion on the drawbacks of this computational 
tool). We hope that future development of numerical realizations of this spe¬ 
cial function in symbolic computational software packages will make its usage 
in the problem under consideration more convenient and reliable. 

Comparing our approach with CKD we conclude that it has both merits and 
drawbacks. On the one hand it is conceptually simpler because it makes use of 
own eigenfunctions for the Smoluchowski operator rather than foreign ones. As 
a consequence to deal with inhnite series of closed expressions for independent 
coefficients is simpler than with inhnite hierarchies of differential-recurrence 
relations where all coefficients are interdependent of one another. Also our 
approach provides the option to deal directly with the probability distribution 
function rather than with statistical moments (i.e., values averaged over it). 
On the other hand it is applicable only to symmetric Maier-Saupe potential. 
Besides, CKD provides a user of the theory with compact analytical formulas 
for MFPT while our approach leads to the expressions that can be treated 
only numerically. 

The approach developed in the present article can not be directly applied to 
the asymmetric potentials, e.g., to that 



ln(n«) 



iQCoffey et al. approach for InCAsr)] 
[t^Coffey et al. approach for ln(A3T)| 


n2 present paper approach 


-0- n3 present paper approach 


1)4 present paper approach 


present paper approach 


present paper approach 


Fig. 4. The dependence of higher eigenvalues 14^ {TXodd in the approach of Cof¬ 
fey, Kalmykov, Dejardin and their coauthors) on the parameter cr for the ordinary 
Maier-Saupe potential (the particular case c = 0 of the general potential (6)). The 
parameter cr defines the barrier height of the potential. 


ln(A„) 




Fig. 5. The dependence of higher eigenvalues A^, on the parameter a on the pa¬ 
rameter a for the ordinary Maier-Saupe potential (the particular case c = 0 of the 
general potential (6)). The parameter a defines the barrier height of the potential. 
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taking place for a uniaxial crystal subjected to a dc bias field. The latter is 
considered in details by CKD in [45], [27], [25]. The reason of failure is purely 
mathematical: odd powers of x in the potential (we remind that cos ip = x) 
hinder the crucial substitution y = x^ (11) leading to the terms of the type ^/y 
that make the resulting equation to be intractable. However there are some 
prospects to generalize the present approach to the potential of the type 


ksT 


(Tsin^ Ip — p,F{cosip, siu'^) — 2cIn sin 


where F(cos'^, siu'^) is an arbitrary function and is the asymmetry pa¬ 
rameter. In what follows we denote F{x) = F (^x, ^1 — x'^^ and F^ means 
the derivative of this function over x. The particular case mentioned above is 
F{x) = x; c = 0. In this case we have SE 


dt 




(1 - X-) fFL 


where Ls is the Smoluchowski operator. The solution for the unperturbed case 
.^ = 0 obtained in the present article yields its eigenfunctions 



A. 


r 


■Vn'’ 


Lsn 


( 2 ) 

m 



V 


( 2 ) 


They are eigenfunctions of the boundary problem and thus form a complete 
system of orthogonal functions so that any function can be decomposed in a 
series over them. In particular we have 

= Y.PkF + T.Qr’tP 

k I 

=j:pTF+j:srF 

i 3 


— \ 

dx - 

— \ 

dx - 


(l-x^) 
(l - x^) 


„(1) p' 
In X 

„( 2 ) p' 
Im X 


Here P^, Q”, P™ and P]" are well dehned coefficients, e.g. 


TDn _ 

— 


dx (^1 — 


,9Vk 


( 1 ) 


-1 


dx 


etc. We seek the solution of SE in the form 

n m 
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After straightforward calculations with taking into account the redesignation 
of dumb indices we obtain the system of equations for the coefficients 




This system resembles the inhnite hierarchies from CKD in the sense that each 
coefficient depends on all others. It can be solved by a perturbation method 
providing that ^ is sufficiently small (.^ << 1). However the present article 
is devoted to the stringent approach to the solution of SE and considering 
perturbative techniques is beyond its scope. 

The general case of the potential H('^, ip) can not be treated by the approach 
developed in the present article. The latter works only if the motion over the 
polar angle ip is separated somehow from that over the azimuthal one (p. 

In conclusion, we show that the approach based on solving the Smoluchowski 
equation via the confluent Heun’s function has a wide range of applicability 
to symmetric uniaxial Maier-Saupe potential from the ordinary case to the 
double-well one with variable barrier width. The solution is uniformly valid for 
any barrier height. Our results for both mean Erst passage time and eigenvalues 
are in good agreement with those of Coffey, Kalmykov, Dejardin and their 
coauthors in the whole range of barrier heights. 
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